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We study the dependence of the CNO nuclear reaction rates on temperature, in the 
^ \ range of 10'' ^ 10^ K, the typical range of temperature evolution from a Sun-like 

star towards a white dwarf. We show that the temperature dependence of the CNO 
nuclear reaction rates is strongly affected by the presence of non-extensive statistical 
\ effects in the dense stellar core. A very small deviation from the Maxwell-Boltzmann 

particle distribution implies a relevant enhancement of the CNO reaction rate and 
' could explain the presence of heavier elements (e.g. Fe, Mg) in the final composition 

of a white dwarf core. Such a behavior is consistent with the recent experimental 
upper limit to the fraction of energy that the Sun produces via the CNO fusion 
■ cycle. 

> 

^ ' 1 Introduction 

M'. 

The classical theory of astrophysical reaction rates is based on the classical 
Boltzmann-Gibbs statistics [1,2]. But, if we check for the applicability condi- 
tions inside a typical stellar plasma (e.g. inside the Sun or white dwarfs) one 
qualitatively obtains that [3]: i) the collision time and the mean collisional 
time are of the same order; ii) the interaction is non-local and the random 
walks are not Markovian; iii) Boltzmann's hypothesis (Stosszahlansatz) is not 
rigorously satisfied; iv) the potential energy is almost equal to the thermal 
energy, i.e. the plasma parameter is not much smaller than one or of order 
unity [4]. 

It can be shown that a more suitable physical description relies on the non- 
extensive statistical mechanics [5,6]. In the following, we will calculate non- 
extensive corrections to nuclear reactions rates of stellar plasma and to their 
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chemical composition. During the evolution of a star from a Sun-like stage to 
the white dwarf stage, at high and very high temperature, the main nuclear 
burning cycle involves nuclei of intermediate weight (like carbon, nitrogen and 
oxygen). Therefore, we are primarily interested in studying the CNO cycle of 
reactions. 



2 The CNO reaction rates and luminosity 



The two most important nuclear cycles in stars, during their evolution, are the 
pp chain and the CNO cycle, whose leading reactions are p-\-p^d-\-e'^ + i'e 
and ^^A^ + p -^^^ O + 7, respectively. The previous reactions account for the 
star's total energy production and luminosity [1]. 

Given a reaction between two reacting nuclei i and j, the approximate Max- 
weUian (M) reaction rate, , is [2] 

M _ 2^/2 ,.-1/2 S r /-, X 

^ij Y (k^T)V2 ^ ® ' 



where and Nj are the particle densities, /ijj is the reduced mass, S is the 
astrophysical nuclear factor and r = 'iEQ/ik^T), Eq being the most effective 
stellar energy (i.e. the energy at which the reactions are most likely to occur). 

If we modify Eq. 1 with the non-extensive energy distribution function, the 
modified non-MaxweUian (NM) reaction rate, rfj^ , can be written as [7] 

'■r='-f(i+f''-5*k|=)-p(-A«). (2) 
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where the "new" most effective energy Eq is found to be the root of the 
equation Eq = £'o[(l + 25Eq/ {}i^T)l'^/^ , and 5 = {l-q)/2 is equivalent to the 
entropic index q. 

The deformed reaction rates rf-^"^ of Eq. 2 give corrections to: i) the luminosity 
fraction due to the pp chain and to the CNO cycle (Section 2); ii) the chemical 
composition of the stellar plasma (Section 3). 
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Figure 1. Log-linear plot of dimensionless luminosity over temperature, for the pp 
chain and the CNO cycle. (Dashed line, S = +0.0045, q = 0.991; dash-dotted line, 
S = —0.0045, q = 1.009). The vertical line shows the Sun's temperature. All curves 
are normalized with respect to the pp luminosity inside the Sun. 

In the following, we choose the value \5\ = 0.0045 for the non-extensive defor- 
mation parameter, in order to work out a few reliable calculations. In Section 4 
we will sketch an argument, based on recent results from solar neutrino ex- 
periments, that we have used to find out the previous value of 5. 

In Fig. 1, the numerical curves of luminosity over temperature are plotted, 
along with non-extensive corrections. It should be observed that: i) the lumi- 
nosity yield of the pp chain is slightly affected by the deformed statistics, with 
respect to the luminosity yield of the CNO cycle; ii) the non-cxtcnsivc CNO 
correction ranges from 37% to more than 53%; iii) above T ^ 2 ■ 10^ K, the 
luminosity is mainly due to the CNO cycle only, thus confirming that CNO 
always plays a crucial role in the stellar evolution, when the star grows hotter 
toward the white dwarf stage. 



3 Chemical composition 



The particle density of each nuclear species involved in the CNO cycle of 
reactions depends on the characteristic time constants of the nuclear reactions, 
that can be defined as 
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Figure 2. Log-linear plot of dimensionless equilibrium concentrations of CNO nuclei 
over temperature. Classical statistics has been used. All curves are normalized with 
respect to the initial density (^'^N)q inside the Sun (see Ref. [2]). 

where H is the hydrogen density and (av) is the thermal average. To first 
approximation, the differential system of equations that describes the CNO 
cycle can be linearized, and it is cast in the straightforward formula 

f = AX, (5) 



where X = {^^C ^^C ^'^N ^^O ^''Of is the density vector, and A is a 5 x 5 real 
matrix whose elements are the time constants r defined in Eq. 4. 

At time i = 0, we choose the characteristic ratios [2] 

C^Oo : C^iV)o : C^O)o = 5.5 : 1 : 9.6 , (6) 



while the densities of all other chemical elements involved in the CNO cycle 
are supposed to be zero. The temperature dependence of a star's equilibrium 
chemical composition, i.e. the solution of Eqs. 5 and 6, is shown in Fig. 2, 
considering classical statistical mechanics. 

If we now perform the previous calculations in the non-extensive statistics 
picture, each nuclear density might be corrected. For the sake of discussion, 
we have plotted in Fig. 3 the corrective factor to the -^^C (that is far the most 

important medium-weighted nuclide) concentration as a function of the plasma 
temperature: the corrective factor is defined as the ratio {}''^C)]sm / i}''^C)M-i 
where (^^C)jvm is the ^^C particle density with non-extensive statistics, while 
(^^C')m is the classically calculated density, always with \5\ = 0.0045. 
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Figure 3. Corrections to the concentration of over temperature, due to deformed 
statistics. (Dashed Une, S = +0.0045, q = 0.991; dash-dotted hne, S = -0.0045, 
q = 1.009). 

We find that the deformed statistics affects the chemical composition of the 
plasma, and that the corrective factor is of order 10%, in the case of ^^C and 
^^O. On the contrary, if we consider the ^^N concentration, the corrective 
factor is about one order of magnitude smaller (and it can be neglected with 
respect to correction of ^^C and ^^O). 



4 Neutrino constraint 



Recent numerical calculations, performed on experimental data from the Sun 
collected over more than thirty years [8], give an upper limit to the ratio 
between CNO luminosity and pp luminosity in the Sun, 

(^) " ^-^^^ ' 



from which the maximum allowable CNO neutrino fiux is, classically, [9,10] 
{<^rNo)M ^ 2.49 • lO^cm-^s-^ . (8) 



The maximum non-extensive neutrino fiux, {^c^o)nm, can be calculated 
through the first order formula [3,11] 

{<^cno)nm ^ {^cno)m exp{-S(3) , (9) 
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where (3 ~ 338.5. Prom Eqs. 7, 8 and 9, one can obtain a suitable value for 
the deformation parameter, namely, 

\5\ ~ 0.0045 or 0.991 < g < 1.009 . (10) 

It should be observed that \5\ ~ 0.0045 is the maximum allowed value for the 
deformation parameter when considering the Sun's interior only. Therefore, 
the maximum value might be still greater when considering extreme physical 
conditions (e.g. high density plasmas in white dwarfs). This value agrees to 
the one derived from kinetic equations in presence of a random microfield 
distribution and given by the expression 5 — 12 F^, where a is a parameter 
related to ion-ion correlation and to the enforced Coulomb cross section [6] . 



5 Conclusions 

We have studied the temperature dependence of the CNO nuclear cycle in the 
range 10^-i- 10^ K. In this range of temperature, the CNO cycle plays a crucial 
role in the stellar evolution from a Sun-like star towards a white dwarf. We 
have found that a small deviation from the Maxwell-Boltzmann energy dis- 
tribution function inside the stellar core implies relevant modifications on the 
nuclear reaction rates, luminosity and chemical composition of the plasma. 
For 5 < (i.e. g > 1), we obtain a remarkable increase of the CNO reaction 
rates with a more relevant contribution to the star luminosity with respect to 
the one obtained in the classical picture. Such a modification is consistent with 
the recent solar neutrino constraints that fix the deformation parameter to the 
value \5\ < 0.0045. The modified CNO reaction rates thus imply a faster evo- 
lution of stellar nuclear plasma towards heavier elements (like Mg and Fc) at 
high temperature. Such a behavior can be very relevant to understand the na- 
ture of matter in white dwarfs stars. In fact, a comparison of current observed 
mass-radius determinations with the theoretical curves seems to confirm that 
the composition of most white dwarfs is dominated by medium weight el- 
ements (carbon and oxygen). However, a small minority of white dwarfs do 
have relatively small radii indicating the presence of iron cores, which presents 
an intriguing puzzle from the point of view of stellar evolution [12]. A slight 
non-extensivity of the system could explain the presence of heavy elements in 
the final composition of white dwarf's core. 
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